Elliptic flow from colour strings 
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Abstract It is shown that the elliptic flow can be successfully described in the colour string 
picture with fusion and percolation provided anistropy of particle emission from the fused string 
is taken into account. Two possible sources of this anisotropy are considered, propagation of the 
string in the transverse plane and quenching of produced particles in the strong colour field of the 
string. Calculations show that the second source gives an overwhelming contribution to the flow at 
^ ' accessible energies. 
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1 Introduction 

The observed elliptic flow in heavy-ion collisions can be conveniently understood if the distribution 
of the observed particles depends not only on the physical conditions realized locally at their 
ly-) ■ production point but also on the global geometry of the event. In a relativistic local theory this non- 

local information can only emerge as a collective effect, requiring interaction between the relevant 
degrees of freedom localized at different points in the collision region. In this sense anisotropic 
flow is a particularly unambiguous and strong manifestation of collective dynamics in heavy-ion 
collisions. The large elliptic flow coefficient V2 can be qualitatively understood as follows. In 
a high-energy collision spectator nucleons are fast enough to move away leaving behind at mid- 
rapidity an almond shaped azimuthally asymmetric overlap region filled with the QCD matter. This 
spatial asymmetry implies unequal pressure gradients in the transverse plane, with a larger density 
gradient perpendicular to the reaction plane (in- plane) . As a consequence of the subsequent multiple 
interactions between the degrees of freedom involved, the spatial asymmetry leads to anisotropy 
in the momentum space [1]. The final particle transverse momenta are more likely to lie in the 
in-plane than in the perpendicular direction, which leads to V2 > as predicted in [2]. 

This general idea has been realized in various mechanisms for the source of elliptic flow. A 
convenient and successful way to describe the flow anisotropy is achieved in the hydrodynamical 
approach, taking into account the unsymmetric shape of the nuclei overlap in the transverse plane at 
values of the impact parameter b different from zero. This description of course assumes collective 
effects to be responsible for the flow. The flow can also be explained in different terms as following 
from the quenching of initially produced particles in the nuclear matter. Since the path inside 
the nuclear overlap of particles emitted along b and orthogonal to it is different, the observed 
dependence on the azimuthal angle 4> of emission is natural. Finally in [3] the authors pointed 
out that the origin of the flow may be already contained in the asymmetry of the initial emission, 
before any collective effects have taken place. They used the Regge picture of particle production 
and showed that the flow may follow from the dependence of the effective emission vertex on the 
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angle between the directions of the emitted particle momentum p and pomeron propagation e. In 
the end the strength of the flow, characterized by the standard coefficient 1)2 was related to the 
model dependent coefficient e in the dependence on (pe) 2 . 

In this note we want to draw attention that the elliptic flow can also be naturally explained in 
the colour string approach, duly generalized to include string fusion and percolation [4, 5, 6, 7]. This 
approach has proved to be quite successful in the description of particle production and correlations 
in the soft part of spectra [8, 9, 10]. In the original formulation in which strings were assumed 
to be just points in the transverse space there was no place for elliptic flow whatever the string 
distribution in the nuclei overlap. However such a picture is of course too simplified. It neglects at 
least two circumstances. First, strings can propagate in the transverse space and so be characterized 
by a vector in this space rather than a point. In the Regge language this corresponds to assuming 
that the pomeron slope a' is substantially different from zero, in contrast to the original assumption 
that the slope is quite small. In the 3-dimensional space the string becomes not orthogonal to the 
transverse plane and one may expect an anisotropic emission of particles in this plane and a non- 
zero elliptic flow. This mechanism is quite similar to the one considered in [3] with the anisotropy 
created already in the initial emission. Of course the actual form of this anisotropy cannot be 
established on purely theoretical grounds but has to be taken in a semiphenomenological manner. 
We study possible forms of this anisotropy and compare it with the experimental situation. Our 
conclusion is that with any choice of anisotropy in the emission of the string the final elliptic flow is 
far too small as compared to the experimental data. So our conclusion is that taking into account 
only the flow from the initial emission, one cannot explain the data. Inclusion of collective effects 
in the form of consecutive interactions seems to be anavoidable. 

This motivated us to study an alternative mechanism of azimuthal anisotropy, using the second 
circumsance characteristic for a realistic string picture, namely that the string is not point-like in 
the transverse plane but occupies a certain non-zero area in the transverse space and that a fused 
string occupies a greater area than the simple string. Then partons emitted at some point inside 
the string have to pass a certain length before they appear outside and are observed. It is natural 
to assume that as they pass through the strong colour field inside the string they emit gluons and 
their energy diminishes. As a result, the particle observed with transverse momentum p outside 
the string has to be born with a greater momentum inside the string, whose value depends on 
the path lengh travelled inside the string and so different for different direction of emission. In its 
spirit, this mechanism is similar to that of particle quenching in the nuclear matter giving the string 
explanation and characteristics of this quenching. Our calculations confirm that a sizable elliptic 
flow follows from this mechanism. Its centrality and transverse momentum dependence agrees with 
the behaviour of the RHIC data [13, 14, 15, 16]. These results also confirm the ones obtained in a 
similar framework using different simplified methods [17, 18] 



2 A string stretched between points (5\ and 

We consider a model in which strings can be formed between different points in the transverse plane, 
say, /?i in the projectile nucleus and /?2 in the target one. Here j3\^ are in fact two-dimensional 
vectors, which we do not specifically mark in the following. The probability to form such a string 
will be given by the distribution in f3 = f3\ — $2 of the soft pomeron 

P{y^) = ^r y e^l^y\ (1) 



2 



y=y' 



Figure 1: The string in the 3-dimensional space. The shown angle is 9 = tt/2 — 9' 

Obviously the pomeron extends to larger distances as y grows, the average distance being < j3 >= 
\/4a'y. This will be also the average length of the string in the transverse plane. 

We shall be interested in the distribution of particles emitted from such a string. To see it we 
have to recall that the string has also some dimension along the z-axis. We assume that the string 
is located in the xz-plane and its direction defines the z' axis in the primed coordinate system 
(V, y',z'). The y' axis is assumed to coincide with the y axis in the original system and the x' axis 
is taken to be orthogonal to y' and z' axes (see Fig.l). 

Our initial assumption, borrowed from the Schwinger mechanism of particle emission in the 
external field, is that particles are emitted isotropically only in the plane orthogonal to the string 
direction. So in the primed system we find that the probabilty to emit a particle with momentum 
p' is 

fi'(y,p') = 6(p' z )L-W* + P' 2 y\ (2) 
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To find the distribution in the original system we have to express p' via the momentum p in the 
original system, that is to make a rotation by angle 9' between the two systems. We find the 
distribution of emitted particles in p± in the original system as 



Angle 0' is given by 



H(V,P±) = -^-e-^ 2 y + ^ cos2e '\ (3) 

TT COS V 

9' = arctan-^-. (4) 

As we have seen, on the average /3 ~ y/a'y and grows with energy of the string W = e y l 2 . The 
length in z can be estimated as 1/W = e~ y / 2 . So angle 6' is quite close to ir/2. Putting 9' = | — 9 
we find that 9 is small and diminishes with energy: 

„ Az e ~ v/2 

9 --fT~-n=- 5 

P \>a'y 

This means that unless a' = the string only radiates in the direction of y, so that 

Kv,P±) = \jls(p x )e-^ 2 y = ^ls(p x )e-^-V-\ (6) 

On the other hand, if a' = then 9' = 0, 9 = ir/2 and we return to the old case of string 
characterized by a single impact parameter with the emission probability 

»(y,P±)*'=o = -e-^. (7) 

TT 

Obviously with a' ^ emission is anisotropic and we expect a non-zero elliptic flow effect. 

One has to understand that this ideal case may have little to do with the realistic behaviour of 
the string, which does not exactly conform to the Schwinger mechanism, since its dimensions are 
finite. To come somehow closer to reality we may soften our assumption that the string does not 
emit particles along its axis at all. Instead we may assume that emission along the axis is simply 
different from emission in the transverse plane, so that the three-dimensional distribution in the 
primed system is 

^(y,P / ) = v / ^( i ) 3/2 ^ (,?p1+p " +p " ) > (8) 



where i] shows the difference in emission along the string and orthogonal to it. If 77 > 1 the emission 
along the string is damped. In the limit r\ — > 00 we return to the previous case of an ideal string. If 
77 < 1 then, on the contrary, emission along the string direction is enhanced. The same distribution 
in the original system is easily obtained expressing the primed momenta via the ones in the original 
system. The final distribution in the transverse momentum p± is obtained after integration over 
p z and we find 

fi(y,P±) = J dp z p(y,p) = ^-^ ex P ( - £(P? + 7Px))> ( 9 ) 
where the final excentricity parameter 7 is 

^ cos 2 9 + £ sin 2 9 ^ ^ 
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and we used 9' = it/2 — 9 

In the high-energy limit when fl->Owe now find 7 = 17 and 



(J-(y,P±) = -y/rjexp \ - dp^ 



If a' = we have = tt/2, 7=1 and 



v(y,p±) = -exp 

7T 



Obviously with 7 different from unity the anisotropy is fully determined by its value. 



(11) 



(12) 



3 Elliptic flow from independent strings 

Once emission from a string is anisotropic, it is obvious that we shall find a non-zero elliptic flow 
in AB scattering. The derivation closely follows [3]. Let the probability to find a string attached 
to point /3i in the projectile nucleus A to be 



Pa(Pi 



d 2 biT A (biWi - h), 



(13) 



where Ta is the profile function, and a similar probability Pb{P2) to find a string attached to point 
Pi in the target nucleus 

Pb{P2) = J d 2 biT B {b - b 2 )t{pi - b 2 ), (14) 

where b is the impact parameter of the collision. The probability to find a string stretched between 
the points Pi and Pi will be determined by the factor p{Pi — Pi), where 



P(P) 



1 



■ exp 



(--£-)■ 



Aixa'y 

So the final distribution of emiitted particles will be 



M(b,p) = J d 2 Pid 2 p 2 P A (Pi)P B (Pi) P (Pi - Pi)fi(y,P,p±), 



(15) 



(16) 



where \i is given by (9) with the cc-axis along the direction of the string, that is, along P = Pi — Pi 
One can transform (16) to the form introduced in [3]: 



where obviously 

and 
with 



M(b,p)= J d 2 pT int (b-P)T(P,y, P± ), 

T int (b) = J d 2 PiP A (Pi)P B (P-Pi) 
T{P,y,p±) =p(P)n(P,y,P±), 



v(P,y,p±) = -^exp 



i_(pPY jppy 

02 +7 ^ 



P 2 



(17) 

(18) 
(19) 
(20) 
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Taking the final x-axis along the direction of p± we get the expression which can be compared 
with [3] 

F(P, P ±) = ^Pr ex P [ - Tl- ~ ^P 2 ( sin 2 X + 1 cos 2 X )l , (21) 

where x 1S the azimuthal angle between (3 and 

As compared to [3] the role of p 2 and f3 2 are reversed in our approach. In our case the distribution 
T is essentially 

T ~ e -/ 32 //3o-?P 2 ( sin2 X+7cos 2 X ) ) (22) 

where /3g = 4a' y and excentricity 7 is given by Eq. (10) (it depends (weakly) on the string 
dimension /3). In [3] one has 

r ~ e -?P 2 -(/3 2 /A))(sin 2 X+7cos 2 x) ) (23) 

where /3q and 7 are expressed via parameters of the emission vertex in the Regge approach 

«• = *.. ^=^4^' (24) 

with r q ~ ro ~ 1 fm. In [3] 7 depends on p 2 (not weakly). 

However it turns out that after integration over <p both forms are essentially equivalent. The 
main difference in the results comes from the choice of parameters. In particular with the standard 
choice of the pomeron slope a' = 0.2 GeV~ 2 our /3 turns out to be quite small at accessible 
rapidities as compared to the value assumed in [3] and becomes of the same magnitude only at 
y ~ 100 

4 Strings homogeneosly distributed 

4.1 Analytic expressions 

To obtain concrete results we have to specify the distributions of strings inside the nuclei. For 
simplicity we assume strings to be homogeneously distributed in the transverse plane inside the 
nucleus: 

PaW) = *. e{R ^\ (25) 

where N s is the number of strings and Ra is the nucleus radius. We also consider a simple case of 
a collison of two identical nuclei, so that, up to coefficient, Tj nt is just the overlap area of the nuclei 
at distance f3 between their centers: 

T mt (/3) = ^(2C-sin(2C), C = arccos^-. (26) 

For the distribution in angle <p between p and b the coefficiens in front of the inclusive cross- 
section are unimportant. So for our purpose we find 

M(b,p)= f m f 2 * dx(2(-sm(2()e~P 2 /fte-tP 2 ( s ™ 2 x+~f cos2 x) i (27) 
Jo Jo 

where 

p m = 2R A + b, P = 4a'y, C = arcc OS t^, ( 28 ) 

ZKa 
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7 is given by Eq. (10). and £ is the inverse string tension which determines the distribution in 
momenta. The dependence of M(b,p) on comes from function (. If we put x = x' + 4> then 



Jb 2 + B 2 - 2bB cos y , . 

C = arccos — — 29) 

2R A v ; 



and we get 

r/3 m 



M(b,p) = f d X '(2( - S m(2C))e- /32 / /3 oe- a P 2 ( sin2 ^' + ^ + ^ cos2 ^' + ^), (30) 



where integration over \' should be restricted by the condition 



lx l< arccos P — 4 (31) 



unless the argument of the arccos is less than —1. 



4.2 Numerical results 

We take the standard value for a' = 0.2 GeV -2 . For the string tension £ we choose £ = 0.25 GeV~ 2 . 
The only left parameter is excentricity 7. From the start one finds that with values of 7 greater 
than unity we obtain negative values for the elliptic flow parameter V2 and with values of 7 less 
than unity we get positive values of V2- Should the elliptic flow come only from the discussed effect 
of string extension in the transverse plane, the experimental data would exclude 7 > 1 and thus 
the naive picture in which the string only emits particles in the plane transverse to its direction, 
as in the Schwinger picture. Rather we have to admit the opposite: the string emits particles 
predominantly along its direction, which leads to 7 < 1. This does not look too exotic if the length 
of the string in the z-direction is much smaller than in its transverse direction. 

However both with 7 < 1 and 7 > 1 the magnitude of the elliptic flow turns out to be quite 
small at accessible rapidities. We studied Au-Au collisiona at rapidity y = 10. In Fig. 2 and 3 
we show our results with 7 = 0.1 which lead to positive values for V2- In Fig. 2 we show V2 as 
a function of impact parameter b at different values of p. In Fig. 3 we illustrate the p dependence 
of the elliptic flow at b = 6. As we see the qualitative behaviour of both b- and p- dependence of 
the elliptic flow agrees with the experimental findings. However values of t>2 at not too peripheral 
collisions are ten to twenty times lower than the data. Only in the limiting case of a very small 
overlap these values become comparable to the observed ones. 

In (4) we show the elliptic flow for 7 = 10, when the resulting V2 is negative. Its values turn out 
to be practically independent of p. The magnitude of | ^2 1 an d its behaviour with b are very similar 
to what we have obtained for values of 7 smaller than unity. Again \v2\ reaches value of the order 
of several percents only at very peripheral collisions. 



5 String fusion as a source of elliptic flow 
5.1 The model 

As we have seen, propagation of strings in the transverse plane gives a certain contribution to the 
elliptic flow, but it is too small compared to the experimental data and can only be noticeable at 
extremely peripheral collisions. So we pass to another source of elliptic flow related to the process 
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Figure 3: Elliptic flow coefficient V2 with 7 = 0.1 as a function of p at b = 6 fm for Au-Au collisions 
at y = 10. 
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Figure 4: Elliptic flow coefficient \v2\ with 7 = 10 as a function of b for Au-Au collisions at y = 10. 
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of string fusion and mentioned in the Introduction. In fact this process is of the same nature as 
implied in the hydrodynamical approach, a sort of collective interaction. The difference is rather 
quantitative, since fusion of strings brings the nuclear matter in the overlap to the liquid state 
gradually, depending on the so-called percolation parameter p, which is the string density in the 
transverse space. Only at values p > 1.1-=- 1.2 drops of nuclear matter (fused strings) on the nuclear 
scale begin to form and at very high values of p they tend to occupy all the overlap space. 

The elliptic flow may have its origin in the fact that fused strings are not symmetric in the 
transverse plane and therefore may emit particle with different probability in different azimuthal 
directions. As a source of this asymmetry one may consider quenching of emitted particles as they 
pass through the fused string area. In this case if a particle is emitted from the forward surface of 
the fused string it will have a greater momentum than the one emitted from the backward surface, 
the latter having to pass through the string loosing its energy. Generally, if the emitted particle has 
to propagate inside the string a path of length I we may assume the probability to see the particle 
of momentum p outside the string to be proportional to 

P(p,l) = e-&1, (32) 

where a characterizes the loss of the transverse energy per unit length. If the fused string is not 
symmetric, then the length I will depend on the direction of emission, that is on the angle between 
the x -axis related to the string and momentum p. 

By itself this anisotropy of particle emission from fused strings cannot give any elliptic flow. 
If the overlap is azimuthally symmetric (impact parameter 6 = 0) then obviously strings will 
occupy arbitrary directions respective to b and the resulting distribution in momenta will be fully 
isotropic. However if the overlap is not symmetric in the azimuthal angle, then fused strings of 
different directions will form with different probabilities. This will give rise to a non-zero elliptic 
flow. A clear limiting case is when p >> 1 and all strings fuse into one which occupies all the 
overlap area and so is unsymmetric if the latter is. 

Below we shall try to develop a more quantitative way to calculate the elliptic flow from string 
fusion. As a typical form of the fused string we shall assume a symmetric almond simlar to the 
shape of the nuclear overlap for collisions of two identical nuclei with radius R at distance b between 
their centers. It is described by the equation 

(x±^) 2 + y 2 = i? 2 (33) 

where x-axis is directed along its minor axis. The two axes of the almond itself a and d with a > d 

are 

a = V4R 2 - b 2 , d = 2R-b, b<2R. (34) 

We shall be interested in the emission of particles at angle x to the minor axis of the almond 
(x-axis) in the forward direction. Because of symmetry we can assume < x < 7r /2- Particles 
emiited at point (xo,yo) in this direction move along the line 

y = ax + c, c = yo — axo, a = tanx (35) 

(see Fig. 5). 

Points x\ and X2 on the surface of the string from which the particle is emitted in the forward 
and backward directions respectively are obtained as solutions of the system of equations (35) and 
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Figure 5: The almond of the nuclear overlap and the path of a particle emitted from point (xo,yo) 
at a given angle to the minor axis of the diamond 
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(33) with signs '+' and '— respectively. They may lie either on different sides of the almond or on 
the same side depending on the values of c and a Trivial manipulations give 



A 2 

and 
where 



+ y / Q^), o(<)-|, (37) 



A 2 = l + a 2 , Qj = A 2 i? 2 - (c- afe/2) 2 , Q 2 = A 2 i? 2 - (c + a&/2) 2 . (38) 
The length travelled by the emitted particle in the forward direction is given by 

1 = \Jxi- x ) 2 + (yi - ?/o ) = A(xi - x Q ) = ^(yi - y ). (39) 

The emission probability will be given by (32). It will depend on the initial emission point 
(xo,yo) and angle %■ To find the total emission probability at given p and \ one nas to integrate 
over all (xo,yo) m the almond area. Up to a constant factor 

f tp 2 
P(P,X)= dx dy e i-^U^o.fo) . (40) 

The total emission probability at given x will be given after integration over all p. Again up to a 
constant factor 

P(X) = J dx dy (l - al(x, x , yo)) • (41) 

Note that the expression for I considerably simplifies for the extreme angles x = an d X = 7r /2- 
If x = then A = 1 and 

b 



Xl = -x 2 = ^R 2 -y 2 --. (42) 
If x = 7r /2 then A/a = 1 and 

2/1 = -V2 = ^R 2 -(\x \+b/2) 2 . (43) 
In these cases one can obtain P{x) m & n analytical form: 



-2 
-3 

and 



P( X = 0) = n - a[^a(2R 2 + b 2 ) + W2 2 (2C + sin(2C))] (44) 



P( x = tt/2) = n - ^cr[8i? 3 - 66P 2 + ft 3 / 2 ! , (45) 
where 

C = arccos — (46) 

and f2 is the almond area 

tt = P 2 (2C-sin(2C)). (47) 

As mentioned, this anisotropy of emission from an unsymmetric string has to be combined with 

the anisotropic distribution of strings in the anisotropic overlap. Choose the system in which the x 
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axis is directed along the impact parameter. Let the minor axis of the almond corresponding to the 
fused string be directed at azimuthal angle 9. Then particle emitted at angle x respective to the 
minor axis of the almond will be emitted at azimuthal angle <j) = x + 9 respective to the direction 
of the impact parameter. It is the (^-dependence which is measured experimentally, so that in our 
formulas for the emission probability and length I one has to put x — ► X ~ 6- The fused string can 
be generally directed at any angle 9 respective to the direction of b with the probability D(Q) So 
the final distribution in the transverse momentum of particles emitted from thsi fused string will 
be obtained as 

Jj^-cjmimPtp,*-*) (48) 

and the distribution in 4> as: 

^ = CiJ dOD{e)p{<t>-e). (49) 

Obviously if the distribution of strings in 9 is isotropic, then integration over 9 will eliminate any 
dependence on <j). 

If the overlap is fully symmetric (b = 0) then obviously we cannot expect any dependence of 
the string distribution in 9. However if it is not then such a dependence certainly arises. 

To see this, take the percolation parameter p to be very large, so that practically all strings 
in the overlap area fuse into a single cluster, which occupies the whole overlap volume. With an 
uasymmetric overlap the only possibility to form such a cluster is to direct its major and minor 
axes along the mjor and minor axes of the overlap. In other words the distribition of strings in 9 
acquires a <5-function form 

D{9) = 5{9) (50) 

and then 

s^ = CF M- < 51 > 
% - CiF (« < 52 > 

where the almond occupies the whole overlap area, so in P one has to put R = Ra, the radius of 
colliding nuclei, and b the actual impact parameter. The obtained asymmetry is quite similar to 
the one in the well-known mechanims in which the emitted partcles are suppressed by final-state 
interactions as they pass through the overlap area [11]. 

The string fusion mechanism allows to obtain something more: the dependence of the sup- 
pression on the string density through the percolation parameter taking not so high values. Then 
strings occupy only part of the overlap area given by factor 

F(p) = 1 - e~ p . (53) 

This gives them a possibility to form different clusters and have different directions in the overlap 
area. If p and so the string length / become small all directions become equally possible, so that 
the distribution in 9 becomes flat and with that the distribition in (f>. 

Unfortunately to study this p-dependence accurately one has to recur to a very complicated 
numerical Monte-Carlo methods to be able to see forms and directions of the string clusters in 
detail. 

To avoid this we may use a crude simplified approach. We may assume that strings (fused 
and simple) on the average are distributed homogeneously in the overlap area. Again on the 
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average, an emitted particle has to pass through the string matter length / which is smaller than 
the corresponding length at p » 1 by factor 



K 1 (p) = y/F(j>). (54) 

The average string tension will be greater than for the single string by factor n{p)~ l where 



W 1 ^-^ (55) 

The distribution in the transverse momentum will then be given by the same formulas (51) and 
(52) with the appropriate rescaling of I and £ 



da 



pdpd(f> 



CPlp,4>,£ ->• £/e,Z ->■ Kil) = C dx dy e i-^itri't*,^.™) (56) 



^ = CiP((/>,a y clk, I y ml) = J dx dy (l - K 1 (p)al(4>,x ,yo)J (57) 

This approach will be referred to as the 'averaged model'. 

Note that distribution (56) formally may be interpreted as following from the dependence of 
the percolation parameter p on the point inside the overlap. Indeed one can rewrite (56) in the 
form 

= C [ dx dyoe- Kp efftt' xo ' y °K p2 , (58) 
pdpdcp J 

where the point- and direction-dependent effective percolation parameter is determined by the 
obvious relation 



/ Pefffaxo vo) = _Jj L= _ am4> y (59) 

y i _ e -w/(*io,w) y/i _ e - P VFW u ' yu; v ; 

In this way this picture coincides with the one proposed in [17] and may be considered as its 
justification. 

Calculations become especially simple if we neglect higher harmonics in the expansion in cos(rn^) 
and restrict to n = 2. Then we can obtain the elliptic flow coefficient by comparing emission at 
angles </> = and 4> = ir/2. In this case for the distribution integrated over p we get an explicit 
expression 

2V2 _ P(0) - P(*/2) 



P(0 + P(^/2) 
I e(p) [8 - 66 + b 3 /2 - 2a(2 + b 2 ) - 3b(2( + sin(2C) 



2 2 - e(p) [8 - 6b + 6 3 /2 + 2a(2 + b 2 ) + 36(2( + sin(2C)) 



(60) 



where 



€{p) = ~xT (61) 



and we have put Ra — > 1 everywhere else, so that 



n = 2( - sin(2C), C = arccos | (62) 
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and a and b are to be measured in units Ra- 

A more elaborate picture ('elaborate model'), leading to a non-trivial ^-dependence of D(9), 
can be based on the assumption that the macroscopic string cluster has the same almond shape 
as the overlap itself but its dimension is reduced by factor k. Then its direction may vary and the 
distribution of strings in 9 will not be so sharply peaked at 9 = Accordingly we assume the length 
of the major axis of the string is 

a = KiyjAR\ -b 2 = Kia (63) 

Its direction respective to the direction of b will be given by r = 9 + ir/2. Symmetry of the 
almond dictates that the distribution does not change if 9 — > — 9 and 9 — > -it — 9. So it is sufficient 
to consider the case when the major axis of the string has its angle r contained in the interval 
between and ir/2 which corresonds to — ir/2 < 9 < 0. Values of D{9)) at other angles will be 
obtained by symmetry. 

The probability that the string of dimension I and major axis forming angle r with the direction 
of b is proportional to the part of almond area A in which such a string can lie. This area is given 
by the integral over the almond area over points x and x' such that they lie on a line forming angle 
r with the x-axis and with |x — x'| = I. From (39) we find 

A = j J d 2 xd 2 x'5(x' -x- l/X)8{y' -y- al/X) (64) 

where a = tanr and A = y/l + a 2 . The average length I of the string will be given by Eq. (63). 
The two-dimensional integral (64) can be easily calculated numerically. Note, that as follows from 
these numerical calculations, with the growth of b and p, the distribuition proportional to A rapidly 
takes the form of the ^-function, concentrated at r = tt/2, that is at 9 = 0, which returns us to 
the previously considered averaged approximation. So using this more elaborate picture is only 
reasonable at relatively small values of b and p 



5.2 Numerical results 

Our model contains a single new parameter a which characterizes the loss of energy in passing 
through the string field. In fact by dimensional reasons a is proportional to the nuclear radius Ra 
Taking R A = RqA 1 ^ we find 

a = aoA 1 / 3 (65) 

where do is a dimensionless and A independent parameter to be extracted from the experimental 
data. 

Values of the percolation parameter p corresponding to given values of the impact parameter b 
were taken from [12] for Au-Au collisions at 62.4 and 200 GeV. They are shown in Fig. 6. 

Calculations of the elliptic flow coefficient v-i from events integrated over the transverse param- 
eter as a function of impact parameter b according to the averaged formula (57) give the results 
shown in Fig 7 for Au-Au collisions at 62.4 and 200 GeV. Calulations according to the same 
averaged picture (Eq. (56) of the transverse momentum dependence are presented in Fig 8 for 
central, mid-central and peripheral collisions. In both cases quenching parameter do was adjusted 
to the mid-central results for events integrated over transverse momenta. The adjusted value was 
do = 0.09 

Passing to our elaborate model, with a non-trivial distribution in 9, we, as mentioned, have 
found that with the growth of p this distribution rapidly takes the form of a 5 function. This is 
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Figure 8: v<i as a function of transverse momentum for Au-Au collisions at 200 GeV according to 
Eq. (56) for peripheral (uppermost curve), mid-central and central (lowest curve) collisions 
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Figure 9: The distribution D(9) as a function of r = tt/2 — 9 for 6 = i?^ and p = 0.5 (the leftmost 
curve), 1.0 (the middle curve) and p = 2. (the rightmost curve, of which only the lowest part is 
shown) 

illustrated in Fig. 9 where we plot this distribution as a function of angle r = it/ 2 — 9 for b = Ra 
and p = 0.5, 1 and 2. As one can see for p = 2 the distribution is completely concentrated at 9 = 
Turning to the data for b and p in Fig 6 we find that for all their values for Au-Au collisions 
at 200 GeV the distribution D(9) is practically given by 8(9) so that our more elaborate model 
gives nothing new as compared the previous crude one. Calculations for Au-Au collisions at 62.4 
GeV, for which the distribution D(9) is well different from a (5-function give however results which 
are also practically identical with the averaged model. This is illustrated in Fig. 10 in which we 
compare V2 for events integrated over transverse momenta as a function of b for Au-Au collisions 
at 62.4 GeV calculated according to Eqs. (49) and (64) Eq. (57) with eg = 0.09. The difference is 
neglegible. So the sophistication implied in the our elaborate model turns out to have no practical 
value. 

On the purely theoretical level one expects that the elaborate model places more emphasis on 
string fusion, so that it should give much less (j) dependence at low values of p below the percolation 
threshold p ~ 1.1 1.2. This is indeed so as shows Fig. 11 in which we plot V2 for averaged and 
elaborate models at fixed b = Ra as a function of p at compartively low values. One observes that 
at values of p below the percolation the non-trivial model gives V2 substantially smaller than the 
averaged one. 
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Figure 11: «2 calculated in the averaged (solid curve) and elaborate models for Au-Au collisions 
fixed b = Ra quenching parameter ctq = 0.09 and p changing from 0.2 to 1.8. 
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6 Conclusions 



We have demonstrated that the colour string model with fusion and percolation can succesfully 
describe the observed elliptic flow in high-energy heavy ion collisions. An important ingredient 
in this desription is anisotropy of the string emission spectra in the azimuthal direction. This 
may follow both from the string propagation in the transverse palne due to a non-zero pomeron 
slope and from quenching of the emitted partons in the strong colour field inside the string. We 
have found that the first source of anisotropy plays a minor role at accessible energies due to the 
fact that the distance travelled by the string in the transverse palne turns out to be small. The 
second source of anisotropy however gives rise to anisotropy, which, upon adjusting the parameter 
of quenching, allows to describe the data quite well both in their centrality dependence and their 
transverse momentum dependemce. 

Our results have been obtained under some substantial approximations. In the simplest case we 
used the over-all averaged picture both as to the form of the string clusters as to their distribution 
in the nuclei overlap. In the more elaborate case we fixed the geometric form of the leading cluster 
and then found its distribution in the overlap. Both approximations have led to practiclly the same 
results . Still carefull comparison demonstartes that the leading cluster approximation gives less 
flow at small values of the percolation paramter, which indicates that string percolations is the 
most important source of the flow. The results obtained here are similar to perevious calculations 
in [17, 18] in the same framework of percolation of strings under different approximations. 

More accurate calculations of the flow in the string percolation model are only possible in 
the developed Monte-Carlo approach. They do not seem simple since one has to find an overall 
quenching for a given distribution of string clusters in the overlap. We plan to conduct such 
calculations in future. 
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